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Abstract 

The multiple scattering interferences due to the addition of several contigu- 
ous potential units are used to construct composite complex potentials that 
absorb at an arbitrary set of incident momenta or for a broad momentum 
interval. 

PACS: 03.65.Nk, 34.90. +q 

Complex absorbing potentials are an important tool in stationary or time dependent 
scattering calculations [||]. They avoid spurious edge effects of the finite "box" where the 
system is enclosed for numerical purposes in wave packet calculations P]. They have been 
also used in other contexts, such as time-independent approaches to reactive scattering 
0-^, or calculations of the microcanonical cumulative reaction probability, and transition 
probabilities in time dependent fields 0,0. Allcock, and other authors have claimed that 
no perfectly absorbing potential can exist, even for one single incident momentum, with 
finite spatial support p. But counterexamples have been found, i.e., potentials that absorb 
perfectly at one selected momentum for an arbitrarily small support do exist ||^. In actual 
collisions, however, the wave packets may have broad translational momentum distributions 
and be peaked around different momenta. Such a case arises, for example, as a result 
of internal state energy discretization of the collision products [TU|. So the challenge is to 



model complex potentials that absorb perfectly at a discrete set of energies, or with sufficient 
efficiency for broad momentum intervals. 

In this letter we shall make use of the interferences between paths associated with "multi- 
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pie collisions" in composite barriers to construct complex potentials that absorb at a selected 
set of incident momenta or for a broad momentum range. Two different complementary 
methods are described and demonstrated. The first method works by successive addition 
of one potential unit for each absorbed momentum and leads to perfect absorption for the 
selected momenta. In the second method no perfect absorption is achieved but it is nu- 
merically much more robust than the former, and allows for efficient absorption in a broad 
momentum range. 

The first construction method is now described. Assume two complex potential units Vi 
and V2 of contiguous, finite supports, Vi being the first one from the left, and V2 the second. 
Let Tp^ and i = 1,2, be the complex transmission and refiection amplitudes for left 
(/) and right (r) incidence for the i-th (isolated) potential unit, and T^'^ or i?'''' (without 
subscript) the corresponding amplitudes for the compound barrier, V = Vi + V2. These 
last quantities may be obtained in terms of the former by the "multiple collision" technique 
Tl| , i.e., by considering the sum of the amplitudes for all the possible "paths" that lead 



eventually to transmission or refiection. In particular. 
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Now assume that 

T[{h) = R[{k{) = Q , (2) 

where /ci > is a particular (dimensionless) wave number, {k = dip/h, where p is the 
dimensional momentum and di the dimensional width of the first barrier). Inserting (^ into 
(|I|) gives T\ki) = R\ki) = 0. In other words, if Vi is a perfect absorber at ki, the total 
potential V is a perfect absorber for ki too. The objective of adding V2 is to absorb also at 
k2 7^ ki. If V2 is naively constructed in such a way that -R2(^2) = ^2(^2) = 0, according to 
the equations (|l|) the composite barrier does not transmit the new momentum, T\k2) = 0, 
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but in general R\k2) 7^ because of reflection in the first barrier, R\{k2) 7^ 0. Instead, if 
the second potential satisfies 

T^(fc2) = 0, (3) 

""^^"'^ R{{h)Rl{h) - T[{h)Tl{h) ' 
total absorption is achieved at A;2, T\k2) = R\k2) = 0. Eq. (^) is obtained by assuming 
R\k2) = and solving for R2 in (|l|) in terms of quantities that depend only on Vi. Even 
though the two barriers alone have non vanishing reflection amplitudes for k2, when they 
are put together the interference effects exactly cancel the global reflection. 

We shall next describe a way to construct Vi and V2 with the desired partial reflection 
and transmission amplitudes. Dimensionless variables will be used throughout. The dimen- 
sionless position x is obtained by dividing the dimensional position by a reference length 
{di in this case), so that the support of Vi has length one when using the x variable. The 
second potential unit, V2, may have a different length, L2 = ^2/(^1. Dimensionless energies 
(kinetic and potential) are obtained by dividing the corresponding dimensional quantities 
m being the mass of the particle. 

Potentials Vi with vanishing reflection and transmission coefficients for ki can be con- 
structed by means of an inversion procedure similar to the one described in The im- 
portant difference with respect to that work is that now transmission is allowed at other 
wave numbers, so that the wave function and its derivative are continuous at x = 1. This is 
necessary to take advantage of interference effects. The boundary conditions to be satisfied 
by the stationary wave function corresponding to an incident plane wave with wave number 
ki at the two potential edges are 

^i(0) = l, ^i(0) = ^A;i, (5) 
^i(l) =7^1(1) = 0, 

where the prime means "derivative with respect to x". To satisfy these four conditions 
the wavefunction between and 1 is written in terms of a functional form with four free 



parameters. By substituting this expression in the four equations (|^), the four parameters 
are determined, and solving in the Schrodinger equation for Vi{x) one finds 

V, = kl + ^'l/^,. (6) 

The simplest choice for ipi is a polynomial, ipi = J2j=o,3(^j^'' ■ From (|^) the coefficients aj 
are readily obtained, 

ao = 1, ai = iki, 02 = — 3 — 2iki, = 2 + iki . (7) 

To construct the potential V2 between 1 and I + L2, or additional units with support between 
two points X = z and x = 2; + L, it is convenient to define the new variable y = {x — z)/L, 
so that the potential unit goes from y = to y = 1. A new wavenumber is also defined 
as k = Lk (remember that L is dimensionless) . The refiection amplitude r\k) for the new 
variables (A;, y) is related to the one for the original variables (/c, x) by r\k) = B} {k)e~'^'^^^ 
. It is much easier to manipulate the constraints (|^) and (Q) using the new set of variables. 



In particular, for obtaining V2, it is first assumed that the wave function '^2{y) corresponding 
to a plane wave incident from the left with wavenumber k2 = i^2^2 obeys the four conditions 

^2(1/ = 0) = 1 + 4 (8) 
%{y = 0) = tk2{l-rl) 
^,{y=l) = %{y=l)=0, 

with rl{k2) = -R2(^2)e~^*'^^ and -^2(^2) given by (H). By assuming, as before, a polynomial 
form, '^2{y) = Z]j=o,3 ^i?/"') substituting in (^), the coefficients are found to be 

bo = 1 + 4, h = ik2{l-rl), (9) 

b2 = -{3 + 2ik2) ~rl{3-2ik2), bs = {2 + tk) + rl{2 - ik) . (10) 

Solving in the Schrodinger equation for the corresponding potential and rewriting the result 
for the original variables, ip2{x) = ^2(1/), one then finds for V2 between x = 1 and x = I + L2, 

V2{x) = kl + fLl. (11) 



The potential Vi + V2 so constructed is a perfect absorber at ki and ^2, and addition of a 
third unit V3 will not change this property, as we have already discussed for the addition of 
V2. By treating the Vi + V2 potential as a new Vi unit, and the new barrier V3 as V2, the 
inversion method can be repeated to build a potential V3 that absorbs k^. This procedure 
can be continued to construct "perfectly absorbing composite potentials" for an arbitrary 
number of momenta . (With little changes the method is also applicable when an infinite 
barrier is put at the right edge of the last barrier. The only difference is that the condition 
"ip'ix = J2Lj) = need not be imposed so that a quadratic polynomial, rather than a 
cubic one, is enough for the last barrier.) Note that the polynomials and the minimal set 
of conditions discussed here have the advantage of providing simple explicit expressions but 
other functional forms for the wave functions may be used, and further conditions may be 
imposed. 

Figure 1 shows the survival probability S{k) = \R\k)\'^ + |T'(fc)p versus k for potentials 
constructed with two units in this fashion (thick solid and dotted lines). The effective 
absorption width around k2 increases with decreasing L2, an effect reminiscent of the the 
broadening of transmission resonances when the walls of a double barrier are narrowed. In 
Figure 2 three units are used to absorb at three different wavenumbers. The prize to pay for 
the requirement of perfect absorption is that, at least for the potential units we have studied, 
numerical instabilities complicate their practical implementation. Even though the potential 
is explicit, the absorption at an arbitrary momentum has to be computed numerically. The 
numerically calculated S at and around k2 becomes very sensitive to small errors in the 
discretization of the potentials if Ti{k2) is extremely small [|13[ (It appears as a denominator 
in the estimate of the error of S caused by errors in way to avoid this problem 

in practice is to truncate the potentials so that S{ki) is not exactly zero (but sufficiently 
close for practical purposes, say S{ki) = 10"^). This increases Ti(fc2) and makes 5* in the 
proximity of k2 much more stable with respect to slight numerical errors due to the discrete 
representation of the potential. However, the difficulties increase when adding more k points. 
A numerically robust alternative is described next. 
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The second method makes also use of interferences between contiguous units, but in a less 
explicit way than the former approach. Now the functional form chosen for the potential is a 
series of equal length complex square barriers with complex energies {Vj}, j = 1,2, ...,N. 
The real and imaginary values of Vj are optimized with standard subroutines according to 
a flexible criterion: The sum of the survivals for a set of s selected points is minimized [111 



f{V,, Vn; h, K) = S{V,, V^; K) . (12) 

a=l 

Note that s and are not necessarily equal. In a generic application the s points are evenly 
spaced in a given interval in order to absorb arbitrary wave packets within the interval. The 
advantage of this functional form for the potential is that, for a given set of values {V^}, the 
total transmission and reflection coefficients and their gradients with respect to variations of 
the barrier parameters can be obtained exactly by multiplication of 2 x 2 transfer matrices. 
These evaluations are very fast, so that many more parameters can be optimized, two for 
each barrier, than for other functional forms |15,p!6|. 



In Figures 1 and 2 the survival probabilities obtained with the two methods for s = 2 
and s = 3 can be compared in the low (dimensionless) wavenumber region, which is the 
important one to minimize the (dimensional) absorbing potential width in the applications. 
The second method provides with just a few barriers excellent survival curves below 0.001, 
which is sufficient for most practical purposes. Note the improvement of the survival curves 
as the number of barriers increases. This method achieves larger absorption widths than a 
previous systematic approach [^|15| in all studied cases [|l^], with the added advantage of 



a very simple implementation. We have also included in both figures the survival curves 
for one of the most frequently used potentials, V = —irjx^, where rj (real) has been chosen 
to minimize S at the two or three selected points in Figures 1 or 2 respectively. A more 
extensive comparison will appear elsewhere ||16|| . 
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Gobierno Autonomo de Canarias (Spain) (Grant No. PB2/95) and Ministerio de Educacion 
y Cultura (Spain) (PB 93-0578). 

6 



REFERENCES 

[1] N. Balakrishnan, C. Kalyanaraman and N. Sathyamurthy, Phys. Rep. 280, 79 (1997). 
[2] R. Kosloff R and D. Kosloff, J. Comp. Phys. 63, 363 (1986). 
[3] D. Neuhauser and M. Baer, J. Chem. Phys. 92, 3419 (1990) 
[4] D. Neuhauser, J. Chem. Phys. 103, 8513 (1995). 

[5] T. Peng and J. Z. H. Zhang, J. Chem. Phys. 105, 6072 (1996); W. Zhu, T. Peng and J. 
Z. H. Zhang, J. Chem. Phys. 106, 1742 (1997). 

[6] T. Seideman and W. H. Miller, J. Chem. Phys. 96, 4412 (1992). 

[7] U. Peskin and W. H. Miller, J. Chem. Phys. 102, 4084 (1995). 

[8] G. R. AUcock, Ann. Phys. 53, 253 (1969); Ann. Phys. 53, 286 (1969); Ann. Phys. 53, 
311 (1969). 

[9] S. Brouard, D. Madas and J. G. Muga, J. Phys. A 27, L439 (1994). 

[10] R. D. Levine and R. B. Bernstein, Molecular Reaction Dynamics and Chemical Reac- 
tivity, Oxford U. P., 1987. 

[11] J. E. Beam, Am. J. Phys. 38, 1395 (1970). 

[12] M. G. Rozman, P. Reineker and R. Tehver, Phys. Rev. A 49, 3310 (1994). 

[13] For any isolated potential barrier the left and right transmission amplitudes are equal 

when the potential level is equal on both sides. 

[14] For some numerical applications, a barrier with positive imaginary part may be prob- 
lematic. In these cases it is necessary to minimize / with the constraint ImVj < 

[15] D. Macias, S. Brouard, and J. G. Muga, Chem. Phys. Lett. 228, 672 (1994). 

[16] J. P. Palao and J. G. Muga, to appear. 

7 



FIGURE CAPTIONS 

Figure 1. Survival S{k) versus k. The thick sohd and thick dotted hues correspond to 
"perfectly absorbing composite potentials" with L2 = 0.5, and L2 = 1.6 respectively (Li = 1 
in both cases), that absorb ai ki — 1 and k2 — 1.2. (For numerical stabihty the theoretical 
Vi is truncated for values of the real or imaginary parts larger than 10^.) The dotted and 
dashed lines correspond to "square barrier composite potentials" with N = 2 and A*" = 3 
respectively, with the same total length as the potential of the thick solid line, and optimized 
for the same values of A; (s = 2). The solid line with circles corresponds to the potential 
—irjx^, where 77 minimizes the sum of the survivals at the two selected wavenumbers. 

Figure 2. S{k) versus k. The thick solid line corresponds to a "perfectly absorbing 
composite potential" constructed with three units to absorb at ki — 1.94, k2 — 4.84, and 

= 7.75, Li = 1, L2 = .008, and L3 = .024. As in Figure 1 this potential is truncated to 
avoid numerical instability. Also shown is the survival for square barrier composite potentials 
optimized for the same points {s — 3), and with the same total length: N — 1, solid line 
; N = 2, dotted line; A^ = 3, dashed line; A^ = 4, dotted-dashed line. The solid line with 
circles corresponds to the potential —irjx^, where 77 minimizes the sum of the survivals at 
the three selected wavenumbers. 
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